A QCD SUM RULE CALCULATION FOR THE y(4140) NARROW STRUCTURE 
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We use the QCD sum rules to evaluate the mass of a possible scalar mesonic state that couples 
to a molecular D*D* current. We find a mass mu*D* = (4.f 4 ± 0.09) GeV, which is in a excellent 
agreement with the recently observed y(4140) charmonium state. We consider the contributions of 
condensates up to dimension eight, we work at leading order in and we keep terms which are 
linear in the strange quark mass rus. We also consider a molecular D* D* current and we obtain 
mn'D* ~ (4.13 ± 0.10), around 200 MeV above the mass of the y(3930) charmonium state. We 
conclude that it is possible to describe the y(4140) structure as a DJI)^ molecular state. 

PACS numbers: 11.55.Hx, 12.38.Lg , 12.39.-x 

There is growing evidence that at least some of the new charmonium states recently discovery in the 
B-factories are non conventional cc states. Some possible interpretations for these states are mesonic 
molecules, tetraquarks, or/and hybrid mesons. Some of these new mesons have their masses very close 
to the meson-meson threshold like the X(3872) [l| and the Z+(4430) ^2,]. Therefore, a molecular inter- 
pretation for these states seems natural. The most recent aquisiton for this list of peculiar states is the 
narrow structure observed by the CDF Collaboration in the decay i?+ Y{A1A0)K'^ J/4>(t>K'^. The 
mass and width of this structure is M = (4143 ± 2.9 ± 1.2) MeV, T = (11.7t|g ± 3.7) MeV @]. Since 
the y(4140) decays into two l'^{J^'-') = 0^(1^^) vector mesons, it has positive C and G parities. 

There are already some theoretical interpretations for this structure. Its interpretation as a conventional 
cc state is complicated because, as pointed out by the CDF Collaboration 0], it lies well above the 
threshold for open charm decays and, therefore, a cc state with this mass would decay predominantly 
into an open charm pair with a large total width. In ref. the authors interpreted the F(4140) as 
the molecular partner of the charmonium-like state F(3930), which was observed by Belle and BaBar 
collaborations near the J/ipLO threshold [5,]. They concluded that the F(4140) is probably a D*D* 
molecular state with J^^ — 0+^ or 2++. In ref. [6| they have interpreted the F(4140) as an exotic 
hybrid charmonium with J^'-'' = 1 ^ . 

In this work, we use the QCD sum rules (QCDSR) 0, B 0, to study the two-point function based 
on a D*D* current with J^'-' = 0++, to see if the new observed resonance structure, y(4140), 
can be interpreted as such molecular state. In previous calculations, the hidden charm mesons 
X(3872), Z+(4430), r(4260), y(4360), y(4660), Z+(4050) and Z+ (4250) have been studied using 
the QCDSR approach as tetraquark or molecular states [10, 11, 12, |3[3i[3j[i11- In some cases a very 
good agreement with the experimental mass was obtained. 

The starting point for constructing a QCD sum rule to evaluate the mass of a hadronic state, H, is 
the correlator function 

n(g) = ^ J d^x e^^-^0\T[jHix)jl,m\0), (1) 
where the current jj^ creates the states with the quantum numbers of the hadron H. A possible current 
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describing a D*D* molecular state with I^J^'-' = 0+0++ is 



(2) 



where a and b are color indices. 

The QCD sum rule is obtained by evaluating the correlation function in Eq. ([1]) in two ways: in the OPE 
side, we calculate the correlation function at the quark level in terms of quark and gluon fields. We work 
at leading order in in the operators, we consider the contributions from condensates up to dimension 
eight and we keep terms which are linear in the strange quark mass m^. In the phenomenological side, the 
correlation function is calculated by inserting intermediate states for the D*D* molecular scalar state. 
Parametrizing the coupling of the scalar state, H = D*D*, to the current, j, in Eq. ^ in terms of the 
parameter A: 



{0\j\H) = A. 

the phenomenological side of Eq. ^ can be written as 



(3) 



A2 



ds ■ 



(4) 



where the second term in the RHS of Eq.Q denotes higher scalar resonance contributions. 

It is important to notice that there is no one to one correspondence between the current and the state, 
since the current in Eq. Q can be rewritten in terms of sum a over tetraquark type currents, by the 
use of the Fierz transformation. However, the parameter A, appearing in Eq. ([3]), gives a measure of the 
strength of the coupling between the current and the state. 

The correlation function in the OPE side can be written as a dispersion relation: 



n 



OPE 



nOPE 



ds 

4m5 ^ ' 



(^) 



(5) 



where p'^^^{s) is given by the imaginary part of the correlation function: ■np'-^^^{s) = Im[n'^-''^(s)]. 

As usual in the QCD sum rules method, it is assumed that the continuum contribution to the spectral 
density, p'^°^^{s) in Eq. (|3]), vanishes bellow a certain continuum threshold sq. Above this threshold, it is 
given by the result obtained with the OPE. Therefore, one uses the ansatz 17 [ 



OPE 



(s)e(s - So) , 



(6) 



To improve the matching between the two sides of the sum rule, we perfom a Borel transform. After 
transferring the continuum contribution to the OPE side, the sum rules for the scalar meson, considered 
as a scalar D*D* molecule, up to dimension-eight condensates, using factorization hypothesis, can be 
written as: 



ds e-^/^^' pO^^(s) , 



4m'i 



(7) 



where 



pO^^(s) = pP^-'is) + p^'^Hs) + p^^'Hs) + p™"(s) + p^'''>\s) + p™"<^"^>(s) , 



(8) 



with 
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where the integration hmits are given by amin = (1 ^ ^1 — -imll s)/2, amax = (1 + — 4m^/s)/2, 
/3min = am\/{sa — ml), and we have used (sga.Gs) = ■ml{ss). We have neglected the contribution of 
the dimension-six condensate (g^G^), since it is assumed to be suppressed by the loop factor I/IGtt^. We 
also include a part of the dimension-8 condensate contributions, related with the mixed condensate-quark 
condensate contribution: 



mix{ss) 



is) 



racm1{ss)'^ 
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(10) 



It is important to point out that a complete evaluation of the dimension-8 condensate, and higher di- 
mension condensates contributions, require more involved analysis p^ . which is beyond the scope of this 
calculation. 

To extract the mass mD-D- we take the derivative of Eq. ([7]) with respect to l/Af^, and divide the 
result by Eq. Q. 

For a consistent comparison with the results obtained for the other molecular states using the QCDSR 
approach, we have considered here the same values used for the quark masses and condensates as in 
refs. O, ini, m, [H, [11, [H, [H, [H: mc(mc) = (1.23 ± 0.05) GeV, rus = (0.13 ± 0.03) GeV, {qq) = 
-(0.23 ± 0.03)3 Qgv3^ (ss) = 0.8(gg), {sga.Gs) = ml{ss) with mg = 0.8 GeV^ {g^G^) = 0.88 GeV^ 
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FIG. 1: The OPE convergence for the Z)*D* molecule in the region 2.2 < < 3.2 GeV^ for ^ = 4.6 GeV. We 
plot the relative contributions starting with the perturbative contribution (long-dashed line), and each other line 
represents the relative contribution after adding of one extra condensate in the expansion: -|- {ss) (dashed line), 
-I- {g'^G'^) (dotted line), -|- mo(ss) (dot-dashed line), -I- (ss)^ (line with circles), -I- m'^{ss)'^ (Une with squares). 



The Borel window is determined by analysing the OPE convergence and the pole contribution. To 
determine the minimum value of the Borel mass we impose that the contribution of the dimension-8 
condensate should be smaller than 20% of the total contribution. 
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In Fig. [T]we show the contribution of all the terms in the OPE side of the sum rule. From this figure 
we see that for A/^ > 2.3 GeV^ the contribution of the dimension-8 condensate is less than 20% of the 
total contribution. Therefore, we fix the lower value of in the sum rule window as M^^^ = 2.3 GeV^. 

The maximum value of the Borel mass is determined by imposing that the pole contribution must 
be bigger than the continuum contribution. In Table I we show the values of M^^^. In our numerical 
analysis, we will consider the range of values from 2.3 GeV^ until the one allowed by the pole 
dominance criterion given in Table I. 

Table I: Upper limits in the Borel window for the D*D* state obtained from the sum rule for different values of 



( GeV) 


M^,,(GeV2) 


4.4 


2.49 


4.5 


2.68 


4.6 


2.87 


4.7 


3.06 



The continuum threshold is a parameter of the calculation which, in general, is connected to the mass 
of the studied state, H, by the relation sq ~ {niH + 0.5 GeV)^. Therefore, to choose a good range 
to the value of Sq we extract the mass from the sum rule, for a given sq, and accept such value if the 
obtained mass is in the range 0.4 GeV to 0.6 GeV smaller than ^/sq. Using these criteria, we obtain sq 
in the range 4.5 < ^/sq < 4.7 GeV. However, because of the complex spectrum of the exotic states, some 
times lower continuum threshold values are favorable in order to completely eliminete the continuum 
above the resonance state. Therefore, here we will also include the result for the D*D* meson mass for 
= 4.4 GeV. 
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FIG. 2: The D^Dl meson mass as a function of the sum rule parameter {M^) for = 4.4 GeV (dotted line), 
^/so = 4.5 GeV (solid line), = 4.6 GeV (dot-dashed line) and ^/so — 4.7 GeV (dashed line). The crosses 
indicate the upper and lower limits in the Borel region. 

In Fig. [21 we show the D*D* meson mass, for different values of ^/so, in the relevant sum rule window, 
with the upper and lower validity limits indicated. From this figure we see that the results are very stable 
as a function of M^. We see also that for ^/sq = 4.4 GeV, we get a very narrow Borel window, and for 
yfso = 4.3 GeV there is no allowed Borel window. 

Using the Borel window, for each value of sq, to evaluate the mass of the D*D* meson and then 
varying the value of the continumii threshold in the range 4.4 < < 4.7 GeV, we get mo-D' = 

(4.14 ±0.08) GeV. 

Up to now we have kept the values of the quark masses and condensates fixed. To check the dependence 
of our results with these values we fix y'so = 4.55 GeV and vary the other parameters in the ranges: 
mc = (1.23 ± 0.05) GeV, m, = (0.13 ± 0.03) GeV, {qq) = -(0.23 ± 0.03)^ GeV^ ml = (0.8 ± 0.1) GeV^ 
In our calculation we have assumed the factorization hypothesis. However, it is important to check how 
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a violation of the factorization hypothesis would modify our results. For this reason we multiply (ss)^ in 
Eqs. ([9]) and by a factor K and we vary K in the range 0.5 < K <2. 

Table II: Values obtained for mn-o', in the Borel window 2.38 < M'^ < 2.72 GeV"^, when the parameters vary 



in the ranges showed. 



parameter 


mn*D* (GeV) 


mc = (1.23 ±0.05) GeV 


4.15 ±0.08 


m, = (0.13 ±0.03) GeV 


4.14 ±0.02 


{qq} = -(0.23 ±0.03)3 Q^yS 


4.14 ±0.03 


ml = (0.8±0.1)GeV2 


4.15 ±0.07 


0.5 < a: < 2 


4.14 ±0.03 



Taking into account the incertainties given above we finally arrive at 

niD'^D; = (4.14 ±0.09) GeV, (11) 

in an excellent agreement with the mass of the narrow structure y(4140) observed by CDF. 
One can also deduce, from Eq. d?]), the parameter A defined in Eq. ([3]). We get: 

A = (4.22 ± 0.83) X 10"^ GeV^ (12) 

From the above study it is very easy to get results for the D*D* molecular state with J^'^ — 0++. 
For this we only have to take rus — and (ss) = {qq) in Eqs. pU]) . This study was already done 
in ref. [lB| considering 4.5 < ^/sq < 4.7 GeV. Although in the case of the D*D* scalar molecule we get 
a worse Borel convergence than for the D*D* scalar molecule, as can be seen by Fig. [31 there is still a 
good OPE convergence for > 2.5 GeV . 
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FIG. 3: Same as Fig. 1 for the D*D* state for ^ = 4.6 GeV. 

If we allow also for the D*D* molecule values of the continuum threshold in the range 4.4 < ^/sq < 
4.7 GeV we get mwD* = (4.13 ± 0.11) GeV. Therefore, from a QCD sum rule study, the difference 
between the masses of the states that couple with scalar D*Dl and D*D* currents, is consistent with zero. 
The mass obtained with the D*D* scalar current is about 100 MeV above the D*D* (4020) threshold. 
This could be an indication that there is a repulsive interaction between the two D* mesons. Strong 
interactions effects might lead to repulsive interactions that could result in a virtual state above the 
threshold. Therefore, this structure may or may not indicate a resonance. However, considering the 
errors, it is not compatible with the observed y(3940) charmonium-like state. 

In Fig. |4] we show the relative ratio (mo'D* — ttio'D')/itid'D* as a function of the Borel mass for 
= 4.55 GeV. From this figure we can see that the ratio is very stable as a function of and the 
difference between the masses is smaller than 0.5%. Although the ratio is shown for — 4.55 GeV, 
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FIG. 4: The relative ratio between the masses of the scalar molecular states mwu' and mo'D* for -^/so = 
4.55 GeV. 

the result is indiscernible from the one shown in Fig. [4] for other values of the continuum threshold in the 
range 4.4 < ^ < 4.7 GcV. 

This result for the mass difference is completely unexpected since, in general, each strange quark adds 
approximately 100 MeV to the mass of the particle. Therefore, one would naively expect that the mass 
of the D*D* state should be around 200 MeV heavier than the mass of the D*D* state. This was, for 
instance, the result obtained in ref. for the vector molecular states DsqD* and DqD* , where the 
masses obtained were: rrij^^^Q, = (4.42 ± 0.10) GeV and mjj^Q. = (4.27 ± 0.10) GeV . 

For the value of the parameter A we get: 

Xjj,jy, ^ (4.20 ± 0.96) X 10-^ GeY\ (13) 

Therefore, comparing the results in Eqs. (|12p and (I13p we conclude that the currents couple with similar 
strength to the corresponding states, and that both, D*D* and D*D* scalar molecular states have 
masses compatible with the recently observed y(4140) narrow structure. However, since the y(4140) 
was observed in the decay F(4140) — > J/ip(f>, the D^D* assignment is more compatible with its quark 
content. 

In conclusion, we have presented a QCDSR analysis of the two-point function for possible D*D* and 
D*D* molecular states with J^'~'' — 0++. Our findings indicate that the y(4140) narrow structure 
observed by the CDF Collaboration in the decay y(4140)i4r+ J /%p(f)K^ can be very well 

described by using a scalar D*Dl current. Although the authors of ref. [1| interpreted the F(4140) as a 
D*D* molecular scalar state and the F(3930) as a D*D* molecular scalar state, we have obtained similar 
masses for the states that couple with the scalars D*D* and D*D* currents. Therefore, from a QCD sum 
rule point of view, the charmonium-like state F(3930), observed by Belle and BaBar Collaborations, has 
a mass around 200 MeV smaller than the state that couples with a D*D* scalar current and, therefore, 
can not be well described by such a current. 

While this work has been finalized, a similar calculation was presented in ref. [20| . However, the author 
of ref. _^20] arrived to a different conclusion. 
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